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We  have  shown  theoretically  that  efficient  multiple  exciton  generation  (MEG)  by  a  single  photon 
can  be  observed  in  small  nanocrystals  (NCs).  Our  quantum  simulations  that  include  hundreds 
of  thousands  of  exciton  and  multi-exciton  states  demonstrate  that  the  complex  time-dependent 
dynamics  of  these  states  in  a  closed  electronic  system  yields  a  saturated  MEG  effect  on  a  picosecond 
timescale.  Including  phonon  relaxation  confirms  that  efficient  MEG  requires  the  exciton-biexciton 
coupling  time  to  be  faster  than  exciton  relaxation  time. 

PACS  numbers: 


Solar  light  would  be  an  important  source  of  clean  and 
renewable  energy  if  the  efficiency  of  inexpensive  solar 
cells  could  be  increased.  Increased  efficiency  can  be 
achieved  through  carrier  multiplication:  Photo-generated 
carriers,  whose  excess  energy  is  greater  than  the  energy 
gap,  can  create  secondary  electron- hole  pairs  via  impact 
ionization  of  the  filled  band.  Through  this  process,  two 
(or  more)  electron-hole  pairs  are  collected  from  each  pho¬ 
ton  instead  of  just  one.  This  process  is  very  inefficient  in 
bulk  semiconductors,  where  impact  ionization  has  a  very 
low  probability  and  carrier  thermalization,  which  always 
competes  with  impact  ionization,  is  much  faster. 


As  first  suggested  by  Nozik,  impact  ionization  may  ef¬ 
fectively  compete  with  cooling  in  NCs,  due  to  the  en¬ 
hanced  rate  of  inverse  Auger  processes  for  carrier  multi¬ 
plication  and  the  “phonon  bottleneck”  suppression  of  car¬ 
rier  relaxation,  leading  to  efficient  MEG  [1].  Soon  after 
this  publication,  Schallcr  and  Klimov  [§]  observed  ultra¬ 
efficient  MEG  by  a  single  photon  in  PbSe  NCs,  using 
band-edge  transient  absorption  measurements.  Later,  ef¬ 
ficient  MEG  was  observed  by  many  groups  using  different 
techniques  in  NCs  of  many  semiconductors:  PbSe 
Si  [^[jCdSe  0-ll|,  InAs  12,  13|,  and 


PbS  |3|,IY||,  »1  M,  Get  be  [9H11I],  InAs  fia,[13jj,  and  m  car¬ 
bon  nanotubes  15| .  At  the  same  time,  some  groups  were 
not  able  to  observe  MEG  in  CdSe  [Hi  and  InAs  14,  17| 
NCs  and  found  that  the  efficiency  of  MEG  measured  in 
PbSe  NCs  [1]  was  appreciably  smaller  than  that  reported 
earlier.  The  diverse  experimental  data  on  the  MEG  effi¬ 
ciency  are  now  converging  to  more  modest  values  for  PbS 
and  PbSe  NCs  [HI,  but  MEG  in  NCs  has  been  shown  to 
be  significantly  more  efficient  than  impact  ionization  in 
bulk  semiconductors  19|, 


Previous  attempts  at  a  theoretical  understanding  of 
the  enhanced  MEG  provide  only  estimations  of  the 
MEG  efficiency  observed  in  NCs  [§,  [2ll-[26|.  A  self- 
consistent  theory  of  this  phenomena  requires  a  currently 
non-existent  theoretical  description  of  both  the  relax¬ 
ation  mechanisms  for  and  couplings  between  the  highly 
excited  exciton  and  multi-exciton  states  in  NCs  [24|.  To 
explain  the  high  efficiency  of  MEG  in  NCs,  the  coherent 


superposition  model  ,  based  on  the  strong  quasi- resonant 
coupling  between  exciton  and  multi-exciton  states  in  a 
NC,  was  proposed  1 3,  24],  Non-coherent  models  for  effi¬ 
cient  MEG  in  NCs  fill  12.‘jI.[2G|  are  based  on  the  important 
observation  that  the  density  of  biexiton  states  is  signif¬ 
icantly  larger  than  the  density  of  exciton  states  at  the 
same  energy  2l| ,  and  the  density  of  trion  states  is  partic¬ 
ularly  important  for  efficient  MEG  [3] .  The  calculations 
of  MEG  efficiency  in  all  of  these  non-coherent  models  are 
based  on  Fermi’s  Golden  Rule,  which  requires  the  final 
biexciton  state  to  decay  much  faster  than  the  rate  of  the 
exciton-biexciton  transition,  an  assumption  that  has  not 
been  justified  either  experimentally  or  theoretically. 


In  this  letter  we  unify  both  approaches  and  consider  a 
single-photon  excitation  coherently  coupled  with  multi- 
exciton-states  in  a  NC  within  a  full  quantum-state  evo¬ 
lution  approach.  The  time-dependent  dynamics  of  our 
modeled  systems  is  described  using  a  large  multiple- 
configuration  basis  representation  of  the  many-electron 
Hamiltonian,  including  energy  non-conserving  exciton 
and  biexciton  decay  channels.  The  calculations  show  that 
even  in  a  closed,  energy-conserving  electronic  system,  the 
excitation  becomes  predominantly  multi-excitonic  on  a 
picosecond  time  scale.  The  initial  single-photon  excita¬ 
tion  is  dispersed  into  the  dense  multi-exciton  state  space 
of  the  NC. 


We  consider  a  single-photon  excitation  of  a  spheri¬ 
cal  PbSe  NC,  where  efficient  MEG  has  been  reported 
M-  The  energy  spectrum  of  electrons  and  holes  at 
the  band  edges  of  four  equivalent  L  valleys  of  bulk 
PbSe  is  described  by  a  four  band  k  •  p  model  [27 1 . 
In  PbSe  NCs  each  electron,  ,  and  hole,  , 

state  within  the  four-band  effective- mass  model  [28|,  is 
characterized  by  the  spatial  angular  momentum  L  = 
0(5),  1(P),  2 (£>),  3(P), ... ,  the  total  angular  momentum 
j  =  L  ±1/2 ,  which  is  the  sum  of  the  angular  momen¬ 
tum  and  the  spin,  the  projection  of  the  total  momentum 
m  =  ±1/2,  ±3/2, ...  ±  j  ,  and  the  spatial  inversion  parity 
7i .  The  energy  levels  of  the  same  symmetry  were  cal¬ 
culated  using  energy  band  parameters  from  Ref.  @  and 
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labeled  by  a  level  number,  ne,h  .  The  four  envelope 
components,  Fi(r)  ,  (i  =  1,2,  3, 4)  of  the  correspond¬ 
ing  single-particle  wavefunctions  can  be  expressed  in  the 
form: 


Fi(r)  =  C^fl7(r)Yff(d^),  (1) 

where  C^'jn'7T  are  constants,  /"(r)  are  specific  radial 
functions  expressed  in  terms  of  the  spherical  Bessel  func¬ 
tions,  LJ  =  L  or  L  +  1  and  rhf  =  to  ±  1/2  (depending 
upon  the  values  of  i  and  7 r),  and  Y~'li  (0,  (p)  are  the 
spherical  harmonics. 

The  Coulomb  coupling  of  an  electron  or  hole  with 
multi-electron-hole  excitations  at  the  same  energy,  which 
are  disallowed  in  the  bulk  due  to  momentum  conserva¬ 
tion,  plays  an  important  role  for  MEG  in  NCs  [24j|.  We 
generate  a  Hamiltonian  matrix  in  the  basis  of  many- 
electron  states  that  are  Slater  determinants  of  single¬ 
particle  eigenstates.  A  many-electron  state  in  this  ap¬ 
proach  is  a  set  of  occupied  conduction  states  (electrons) 
and  unoccupied  valence  states  (holes)  with  an  arbitrary 
but  consistent  sign  convention  for  the  Slater-determinant 
permutations.  We  include  many-electron  basis  states  up 
to  some  chosen  energy  cut-off  that  we  adjust  until  the 
results  converge  (see  Fig.  [[]).  We  modeled  the  Coulomb 
interaction  using  effective  dielectric  constants  for  the  NC, 
ks  ,  and  surrounding  media,  ng  ,  leading  to  the  direct 
Coulomb  interaction  between  electrons  at  positions  r 
and  r’  in  the  NC,  Vc(r,r')  =  e2/{nAr  —  r*'|  ),  and 
the  surface  mediated  Coulomb  potential  [29j,  Vs  (r,  r')  = 
e2a/(kr’\r  —  a2r'/r'2\) ,  where  a  is  the  NC  radius  and 
kT1  =  2 (ks  —  Kg)/ ks(ks  +  2 Kg)  .  We  also  account  for  the 
electron  and  hole  interaction  with  their  image  potentials 
described  by  Ve(r)  =  T4(r)  =  0.514(r,r)  ■ 

The  diagonal  elements  of  this  Hamiltonian  matrix  are 
taken  to  be  sums  of  the  single-particle-state  energies  (in¬ 
cluding  image  potentials)  and  the  Coulomb  interactions 
among  electrons  and  holes  (including  the  exchange  inter¬ 
actions  from  Fermi  permutations) .  Off-diagonal  elements 
couple  basis  states  that  differ  in  either  two  or  four  sets 
of  single-state  occupations  via  the  Coulomb  interaction. 

Because  the  single-particle  states  in  the  spherical  NCs 
are  represented  in  terms  of  spherical  harmonics,  it  is  most 
efficient  to  express  the  Coulomb  potentials  Vn(r,  r’)  and 
Vs{r,r’)  in  terms  of  spherical  harmonics  [29j] : 
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where  T i(M/,f)=Em1lm(fl'^WM).  r  and 
r'  are  represented  by  {r,  $,</>}  and  in  spher¬ 

ical  coordinates,  r>  =  max(r,  r') ,  and  r<  =  min(r,  r')  . 
Using  these  expressions,  we  have  performed  integrations 
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FIG.  1:  Time  dependent  evolution  of  the  2e —2h P/j2  ex¬ 
citation  created  in  the  2  nm  radius  PbSe  NC  with  ks  =  Kg  =  5 
by  a  single  E 7  =  2.95  eV  photon  calculated  for  the  closed- 
system.  The  upper  panels  (a,b,c)  show  the  occupation  proba- 
blities  of  exciton,  bi-exciton  and  tri-exciton  states  over  0.2  eV 
wide  energy  bins  at  three  times.  Panel  (d)  demonstrates  con¬ 
vergence  of  the  time  dependence  of  the  relative  number  of 
excitons  ( Nx{t ))  —  1  ,  defined  in  Eq.  (j3j) ,  achieved  by  adjust¬ 
ing  the  energy  cut-off.  The  vertical  lines  a,  b  and  c  show  the 
times  corresponding  to  the  three  upper  panels. 


involving  the  angular  components  analytically  by  exploit¬ 
ing  the  orthonormality  of  the  spherical  harmonics. 

The  remaining  radial  parts  of  the  integrals 
of  the  Coulomb  matrix  elements  take  the  form 
fdrfdr'f™(r)f!f{r)fZ(r,)ff,(r')[(r<)l/(r>)l+1} ,  or 
f  dr  f  dr' frn(r)f?{r)f£(r')ff’(r')  (1  /rr’)l+1 ,  where 
/™(r)  are  the  radial  functions  defined  in  Eq.  0  for  the 
envelope  eigenfunctions  that  are  expressed  via  spherical 
Bessel  functions.  In  our  computational  implementation, 
we  approximate  these  Bessel  function  integrals  as  needed 
by  Monte-Carlo  sampling  and  store  the  result  for  future 
reference.  A  sufficient  number  of  samples  was  taken 
to  reach  a  specified  level  of  precision  for  the  matrix 
elements  [30).  At  a  4.8  eV  cut-off  used  to  simulate 
the  2  nm  radius  NCs,  the  basis  contains  hundreds  of 
thousands  of  many-electron  states,  and  the  Hamiltonian 
matrix  contains  tens  of  millions  of  non-zero  elements. 

To  incorporate  the  interaction  involving  a  single  pho¬ 
ton,  we  adopt  the  standard  effective-mass  approxima¬ 
tions  that  the  wavelength  of  the  photon  is  much  larger 
than  the  NC  size  and  that  the  contributions  of  the  Bloch 
functions  are  dominant  in  the  integral  over  the  electron 
momentum  operator  [28| .  In  the  PbSe  NCs  with  symmet¬ 
ric  conduction  and  valence  bands,  we  assume  that  the 
photons  are  coupled  exclusively  to  symmetric  electron- 
hole  pairs,  in  which  both  the  electron  and  the  hole  have 
the  same  quantum  numbers.  As  a  result  the  NC-photon 
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FIG.  2:  Time  dependent  evolution  of  the  relative  num¬ 
ber  of  excitons  ( Nx(t ))  ,  defined  in  Eq.  (0,  created  by  a 
single-photon  excitation  of  various  optical  transitions  in  PbSe 
NCs  with  radii  2  nm  (a)  and  3  nm  (b)  calculated  for  the 
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2D  ,  and  IH  respectively  with  energies  shown  in  units  of 
the  effective  energy  gap  Di(a)  .  Panel  (c)  shows  the  evolu¬ 
tion  of  ( Nx(t ))  created  by  a  single-photon  excitation  of  the 
2eDj/2  —  2hD^2  transition  in  a  2nm  radius  PbSe  NC  calcu¬ 
lated  for  various  combinations  of  the  dielectric  constants  of 


the  nanocrystal  ks  and  the  surrounding  media  K,g  . 


interaction  is  characterized  by  a  single  parameter  cor¬ 
responding  to  the  coupling  energy.  Furthermore,  our 
results  are  obtained  in  the  weak-coupling  single-photon 
limit,  where  the  actual  coupling  strength  is  irrelevant 
provided  that  it  is  sufficiently  small. 

Let  us  begin  from  a  study  of  the  closed  system  evolu¬ 
tion.  Starting  from  an  initial  single-photon-state,  elec¬ 
tronic  states  are  excited  and  evolve  through  Coulombic 
interactions.  In  each  simulated  experiment,  the  photon 
energy  was  chosen  to  match  one  of  the  optically  active 
exciton  states,  because  the  NC  absorption  spectra  are 
not  affected  by  the  coupling  between  excitons  and  multi- 
exciton  states  [24].  We  evolve  the  initial  state  using  the 
Schrodinger  equation,  d\^f(t))/dt  —  —  iH\^(t)) /H  ,  where 
H  is  our  Hamiltonian  that  includes  all  interactions  con¬ 
sidered  above.  Using  a  differental  equation  solver  [3l|  for 
integration,  we  adjusted  the  precision  until  we  obtained 
convergent  results. 

Figure  [T|  shows  the  time-dependent  evolution  of  the 
2eF>i/2  —  2hP\%  excitation  created  by  a  single  pho¬ 
ton  calculated  for  our  closed-system  description.  The 
upper  panels  (a,b,c)  show  the  occupation  probablities 
||  (fc|'F(f))  || 2  of  exciton,  bi-exciton  and  tri-exciton  states 


created  from  the  photon  at  three  different  snapshots  in 
time.  The  initially  created  excitons  are  replaced  by  the 
biexcitons  in  approximately  0.3  ps.  Formation  of  tri- 
excitons  at  the  excitation  energy  considered,  2.95eU  ~ 
3.73£d(a)  where  E\{a)  is  the  effective  energy  gap,  is 
limited,  however,  in  our  single  valley  model,  which  al¬ 
lows  only  two  electron-hole  pairs  to  occupy  the  ground 
state.  Figure  [lid)  shows  the  time  dependence  of  the  ex¬ 
pectation  value  for  the  relative  number  of  excitons: 

(Nx(t))  =  J2nH  ll(fc^(i))l|2/ElK^W)H2>  (3) 

n=l  fceTCn  l>  0 


where  K,n  is  the  set  of  n  -exciton  basis  states.  In  the 
small  NC  with  small  dielectric  constant  considered  in 
Fig.  n  the  strong  Coulomb  coupling  rapidly  creates  a 
bi-exciton  with  90%  probability. 

We  also  investigate  how  ( Nx(t ))  depends  on  the  ex¬ 
citation  energy  and  the  Coulomb  coupling,  which  in 
our  model  is  controlled  by  the  NC  size  and  dielec¬ 
tric  constants.  Figure  Ufa)  compares  the  formation  of 
multi-excitons  from  a  single  photon  excitation  of  dif¬ 
ferent  optical  transitions  in  the  2  nm  radius  NC.  The 


2ePljl  —  2»P£  transition  appears  to  provide  a  thresh¬ 
old  for  efficient  MEG  generation  in  PbSe  NCs:  while 
lower  transitions  saturate  at  around  10  fs  to  a  low  prob¬ 
ability  of  biexciton  generation,  the  multi-exciton  number 
continues  to  increase  when  the  threshold  is  reached.  This 
is  because  the  2 eP  and  2 hP  states  are  the  first  that  can 


,1/2 


be  strongly  coupled  with  trions  of  the  same  energy  [24] . 
Increasing  the  NC  size  [Fig.  [U[b)]  and  the  dielectric  con¬ 
stant  [Fig.  ETc)]  significantly  increases  the  time  of  MEG 
generation.  Both  effects  are  the  direct  result  of  the  re¬ 
duced  Coulomb  coupling. 

Note  that  the  evolution  of  ( Nx(t )}  does  not  show 
strong  oscillations  when  MEG  is  efficient,  even  for  a 
closed  system.  Oscillations  are  seen  only  in  the  strik¬ 
ing  cross-over  behavior  in  Fig.  [2]  (a)  and  (b),  which  is  the 
result  of  the  different  coupling  strengths  and  detunings 
between  the  multi-exciton  states  and  the  exciton  created 
by  the  resonant  light.  The  coherent  oscillations  connect¬ 
ing  the  exciton  and  biexciton  states  are  suppressed  due 
to  the  large  number  of  interacting  biexciton  states.  Once 
the  excitation  evolves  into  the  dense  multi-exciton  state 
space,  it  does  not  appreciably  revive. 

Phonon-assisted  assisted  relaxation  of  carriers,  which 
was  not  considered  above,  can  strongly  affect  the  MEG 
efficiency  [24[.  The  evolution  of  ( Nx(t ))  ,  including 
phonon  decay,  is  described  by 


where  rn  is  the  relaxation  time  of  the  n  -exciton  state 
that  for  simplicity  is  assumed  to  depend  only  on  n .  As 
states  decay,  they  are  removed  from  the  Hilbert  space 
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FIG.  3:  Evolution  of  the  average  number  of  decayed-to-the- 
ground-state  excitons  {Dx(t))  created  by  a  single-photon  ex¬ 
citation  of  the  —  transition  of  a  2  nm  radius 

PbSe  NC  calculated  for  various  combinations  of  thermaliza- 
tion  times.  The  inset  shows  the  distribution  of  inverse  cou¬ 
pling  times  connecting  optically  active  exciton  states  and  bi- 
exciton  states,  which  are  equal  to  the  interaction  Hamiltonian 
matrix  elements,  divided  by  h  .  The  1/ps  scale  in  the  inset 
allows  comparison  of  typical  coupling  times  to  thermalization 
times.  The  gray  region  of  the  coupling  strengths  is  unresolved 
at  our  precision  and  excluded  from  the  distribution. 

and  populate  Dn(t )  ,  denoting  each  “ground”  n -exciton 
state  and  governed  by  Dn  =  J2keic  ll(^l1^)ll2/7h  with 
Dn{  0)  =  0. 

To  find  the  relative  number  of  excitons  follow¬ 
ing  decay,  we  take  the  weighted  average  (Dx)  = 
J2n  nDn/  Drn  ■  We  plot  this  average  in  Fig.  ED  for 
various  rn  lifetimes.  The  introduction  of  the  long  and 
equal  relaxation  times  t\  =  T2  =  73  =  1  ps  has  little  im¬ 
pact  upon  the  time  dependent  evolution:  the  system  re¬ 
mains  almost  closed.  Indeed,  these  times  are  much  longer 
than  an  average  time  of  the  exciton-biexciton  coupling  as 
one  can  see  in  the  inset  of  Fig.  [3]  Reducing  the  biexciton 
and  trion  decay  times  to  T2  =  73  =  100  fs  increases  the 
rate  of  MEG  in  agreement  with  Ref.  [24| .  It  is  interest¬ 
ing  to  note  that  the  saturation  value  for  MEG  depends 
only  upon  t±  .  If  this  time  is  shorter  than  the  exciton- 
biexciton  coupling  time,  ~  300  fs,  MEG  is  significantly 
suppressed. 

Our  calculations  of  the  dynamics  and  efficiency  of 
MEG,  unfortunately,  cannot  be  compared  with  experi¬ 
mental  data  measured  in  PbSe  NCs  because  they  were 
based  on  a  single-valley  model.  The  lack  of  a  precise 
knowledge  of  inter-valley  coupling  does  not  allow  us  to 
take  this  effect  into  account.  The  single- valley  approxi¬ 
mation  should  significantly  underestimate  the  density  of 
multiple-exciton  states.  For  example,  in  real  PbSe  NCs 
one  can  have  eight  electron-hole  pairs  at  the  ground  exci¬ 
ton  state  instead  of  only  the  two  ground  excitons  allowed 


in  our  model.  The  inter-valley  coupling  significantly  in¬ 
creases  the  density  of  multi-exciton  states;  the  k  -exciton 
state  density  is  increased  by  at  least  a  factor  of  42fc  for 
PbSe  with  its  4  valleys.  For  this  reason,  MEG  should  be 
even  more  efficient  in  a  model  that  includes  inter-valley 
coupling. 

In  summary,  our  calculations  unambiguously  demon¬ 
strate  that  highly  efficient  MEG  can  be  observed  in  small 
NCs.  The  effect  is  enhanced  by  a  high  density  of  biex¬ 
citon  states  that  are  strongly  coupled  with  optically  cre¬ 
ated  excitons.  Fast  multi-exciton  thermalization  acceler¬ 
ates  the  formation  of  multi-excitons  in  the  ground  state 
which  should  improve  efficiency  of  extraction  of  electron 
hole  pairs  from  NCs. 
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